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QQ . We analize the statistical mechanics of a long-range antiferromagnetic model defined on a D- 

' dimensional hypercube, both at zero and finite temperatures. The associated Hamiltonian is derived 

from a recently proposed complexity measure of Boolean functions, in the context of neural networks 
learning processes. We show that, depending of the value of D, the system either presents a low 
temperature antiferromagnetic stable phase or the global antiferromagnetic order disappears at 
any temperature. In the last case the ground state is an infinitely degenerated non-glassy one, 
composed by two equal size anti-aligned antiferromagnetic domains. We also present some results 
for the ferromagnetic version of the model. 
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I. INTRODUCTION 



(— I , The success in understanding the collective behavior of conventional magnets lead in the past decade to an increased 
Q ■ interest in systems that exhibit a novel type of ordering as a consequence of "frustration" . In particular, geometrically 
O \ frustrated systems (having as a basic building block triangles of antiferromagnetic bonds) have been the object of 
extensive theoretical studies (see [1] and references therein). It was found that frustration gives rise to spectacular 
and often unexpected behaviors at low temperatures due to the high degeneracy of the ground state. Recently, more 
^ , and more experimental realizations of such geometrically frustrated systems have been achieved. 2D and 3D materials 
fSj ' have been assembled from triangles and tetrahedra such that neighbouring geometrical units share either a common 
T-H \ edge or a common corner [2] . Such new materials have stirred a new debate, whether frustration alone is strong enough 
\l • to destroy magnetic long range order. Spin glass behavior is usually associated to the presence of both disorder and 
frustration. It has been argued that frustration alone can lead to large ground-states degeneracies, accompanied by 
' extensive ground-state entropies, but cannot produce a sufficiently "rough" free energy landscape necessary for the 
■ development of a thermodynamical "glassy" state [4,5]. 
"p-, , In contrast with the previously mentioned systems, non-geometrical frustration may arise as an effect of antifer- 
romagnetic long-ranged interactions (e.g., dipolar. Coulomb, etc.). It is well known that the competition between 
short-range ferromagnetic exchange interactions and long-range antiferromagnetic ones in hypercubic Bravais lattices 
give rise to low temperature long range order associated with lamellar structures [3,6]. Although those systems 
' O [ presents a very rich behaviour concerning both its equilibrium [3]- [6] and non-equilibrium [7,9] properties, in all these 
^ ■ cases the ground state presents a finite degeneracy. Even in the case of a 3D system with Coulomb interactions, where 
some evidence of glassy behaviour has been obtained [9], it is of purely dynamical origin, the ground state still having 
a finite degeneracy. 

In this work we analize the zero and low temperature equilibrium properties of a long range antiferromagnetic spin 
model defined on a Z?-dimensional hypercubic cell. The associated Hamiltonian is related to a recently proposed 
; ^ ■ complexity measure of Boolean functions in the context of neural networks learning processes and it is described in 
I section II. Although the original motivation for studying the statistical mechanics of this model is related to the 
the field of Boolean functions complexity and neural networks [10], we found that it presents some very peculiar 
properties that provide new insights about the relationship between frustration and order-disorder low temperature 
properties. Spin models defined on hypercubic cells have been repeatedly used in the past to study both equilibrium 
and dynamical properties of spin glasses [11]- [13]. In that context, the main interest in hypercubic cell models is 
to analyze how their properties change as the dimension (and therefore the connectivity of the lattice) increases, 
thus resembling fully connected models like the Sherrington-Kirkpatrick one [14], whose general properties are rather 
well understood [15,16]. One of the main questions today concerning spin glasses is whether the physical scenario 
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that emerges in long-range fully connected models holds in systems with short range interactions or not. One step 
further towards the answer of that question is to consider systems with long range interactions that decay with the 
distance between spins. In this work we analize the effects of frustrating long range antiferromagnetic interactions 
in an hypercuibc cell without disorder. The interactions in our model decay with the Hamming distance between 
spins, a natural metric in hypercubic cells, instead of the Euclidean distance. Notice that the Hamming distance in 
the hypercubic cell is equivalent to the chemical distance (i.e., the number of steps required to go from one node to 
another along the shortest path). Hence, the present study may also provide some insigths about the possible effects 
of frustrating interactions we expect in models defined on lattices with complex network topologies, such as small 
world networks (see Rcf. [17] and references therein). 

We show that, for some particular values of the dimension D the frustrating long range interaction gives rise to 
a ground state that is neither disordered (i.e., non-glassy) nor completely ordered but, it is composed of two equal 
size anti-aligned antiferromagnetic domains. The ground state is thus partially ordered. Moreover, it has an infinite 
degeneracy in the thermodynamic limit (although this degeneracy grows logarithmically with the number of spins, 
and hence its entropy per spin goes to zero). 

In section III we analyze the frustrated ground state of the system for the particular set of values of D such that 
it exhibit global frustration. In section IV we analize the finite temperature behaviour of the model using Monte 
Carlo simulations. In section V we present Monte Carlo results for the ferromagnetic version of the model, which are 
of interest in the general context of long range interacting systems. Finally, we present some general conclusions in 
section VI. 



II. THE MODEL 



In Ref. [10] a measure for the complexity of Boolean functions of D Boolean variables was introduced in the context 
of feed-forward neural networks learning processes. The 2^ possible inputs of a Boolean function of this type are 
located at the corners of a D-dimensional hypercubic cell. If we associate a spin variable Si = ±1 to each corner i of 
the cell every possible Boolean function corresponds to a particular spin configuration of this system. The complexity 
measure is based on the concept that complex Boolean functions assign different outputs to similar inputs, where the 
similarity is measured by the Hamming distance between them. Thus, in the magnetic analogy, configurations where 
neighbour spins tend to be anti-aligned correspond to functions of increasing complexity. The complexity measure 
can be then mapped into minus the long range antiferromagnetic Hamiltonian [10] 

j^ = E./N)^.^. (1) 

hi 

where is the Hamming distance between spins i and j in the hypercubic cell and 

^ ' \ otherwise ^ ' 

M{r) being the number of nearest neighbors at distance r of a given spin and equals the binomial coefficient 

^^^) = ^^^'^)=(^^- 

The distribution of Boolean functions with high complexity is then associated with the low energy spectrum of 
Hamiltonian (1) which can be characterized by the low temperature statistical mechanics behaviour of the associated 
spin model. The Hamiltonian (1) can be rewritten as a sum of terms, each of them taking into account the interactions 
between pair of spins at a different Hamming distance: 

r=D/2 

H= E S^^o, (4) 

r—1 [i.j I Hamming{i,j)—r] 

Notice that the number of interaction terms in Hamiltonian (4) depends on D; as D increases a new term is included 
every time D takes an even value. For instance, in the cases D = 2 and 3 only nearest neighbor spins interact; in the 
cases D = 4 and 5 an interaction between next nearest neighbor spins is added, etc.. Hence, the antiferromagnetic 
nature of the interactions introduce local frustration for D > 4, in the sense that not all the couplings inside groups of 
spins closely located can be simultaneously minimized. However, this does not necessarily generate global frustration. 
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in the sense of a non-ordered ground state with high degeneracy, due to the particular form Eq.(2) of the interaction 
coupUng J{r). Global frustration appears only when the number of interaction terms included in the Hamiltonian 
(4) is even, that is, when the dimension I? is a multiple of four D = 41, I — 1,2,3, . . .. In all the other cases the 
ground state is the antiferromagnetic one (i.e., the Neel state with every pair of nearest neighbor spins anti-aligned), 
with degeneracy 2 (corresponding to the inversion of all the spins). To see this let us calculate the energy of the Neel 
state. Every term in Eq.(4) corresponding to a fixed Hamming distance r equals in this case (— 1)''A^ M{r), N = 2^ 
being the total number of spins (notice that the interaction between every pair of spin in Hamiltonian (1) is counted 
twice for simplicity; this corresponds to a renormalization of the energy units by a factor two, and does not affect the 
results). Hence, the energy of the Neel state is given by 

r=D/2 

en = n (-ir (5) 

r=l 

which equals zero when D/2 is even and E^/N = —1 when D/2 is odd. Moreover, in the last case is easy to see 
that a single spin flip increases the energy in a quantity AE = 2. Therefore, the Neel state is stable against small 
perturbations. Monte Carlo simulations presented in section IV provide further evidence that the ground state is the 
Neel one. Moreover, the finite temperature simulations show that there is a finite critical temperature below which 
the system presents long range antiferromagnetic order when D/2 is odd. 

In the next section we analyze the frustrated ground state for the case D = 41, I = 1,2,3, 

III. THE FRUSTRATED GROUND STATE 

From exhaustive enumeration of the energy values for the case D ~ 4, where the Hamiltonian contains two terms, 
it was found that the ground state of the system consist of two equal size antiferromagnetic domains with an energy 
value equal to —0.5. Monte Carlo simulations in higher dimensions {D = 4,8,12,16) have also shown that the 
energy value of the ground state of the system is equal to —0.5 and has the same structure found for the case for 
D = 4. (Long relaxation times and very low cooling were used in different simulations, all of them confirming the 
structure and energy value of the ground state). Energy calculations for states composed by 4 and 8 anti-aligned 
antiferromagnetic domains(D = 8, 16) also yield values higher than —0.5. Moreover, as we will demonstrate next, 
the energy per spin of the state composed by two antiferromagnetic domains equals —0.5 for arbitrary values of D. 
Since the antiferromagnetic state has always zero energy, all these result suggests that the two antiferromagnetic 

domains structure has the minimal energy for arbitrary dimension D = 41, I = 1,2,3, The total magnetization 

and staggered magnetization is zero for the ground state. In figure 1 a representation of the ground state for £) = 4 is 
depicted, where the state of the spins is indicated by filled or empty circles. The antiferromagnetic domains correspond 
to two hypercubes of dimension D — 1 each (see Fig. 1). Each spin in this state has D — 1 nearest neighbor spins 
anti-aligned, which belong to the same domain, and one nearest neighbor spin aligned, which belong to the other 
domain. Once we fixed two nearest neighbor spins to be aligned, the rest of the structure is completely determined. 
Since every spin has D nearest neighbors, the degeneracy of the state is 2 D = 2 log2iV, where the factor 2 results 
from the inversion operation of all the spins. 
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FIG. 1. Ground state of the antiferromagnetic Hamiltonian for D = 4 . Spin states are represented by empty or filled circles 
and the bond between spins indicate that they are first nearest neighbors. 



We now demonstrate that the energy per spin Efi of the state composed of two antiferromagnetic domains equals 
—0.5 for any dimension D = 41, I = 1,2,3, can be written as 

n=D/2 

Ed= ^d{n) (6) 

n=l 

Ed{n) being the interaction energy between an arbitrary chosen reference spin and all its neighbors located at a 

Hamming distance n. 

We denote the two antiferromagnetic domains as regions A and B. Consider now an arbitrary spin that belongs 
to the domain A. Spin Si has D — 1 anti-aligned nearest neighbors (NN) spins in the domain A and only one first 
nearest neighbor spin S2 in the domain B. Spin 5*2 has the same orientation than the spin Si. As the total number 
of NN spins in dimension D is exactly D, the first term of Eq.(6) equals 



which can be expressed as 



Let us now consider the next nearest neighbors (NNN) term n 2 in Eq.(6). Note that the NNN spins of Si in 
the domain B (anti-aligned to Si) are the D — 1 spins NN of ^2 in B; the remaining spins NNN to Si in A have the 
same orientation of it. Since the total number of NNN spins is C{D, 2) we have that 

F (2\ = J-R^ + C{D,2)-{D-1) 
'^^ ' C{D,2) C(D,2) 
_ C{D,1)-C{D,0) 
~ ^ C{D,2) +^ 

Following the same analysis the 3rd term of Eq.(6) equals 
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and the n — th order term can be written as 



p / x r, I]fc=o(~l) C{D,k) 

Equation (11) can be rewritten as 

£<i(n) = (-l)"+i(2a(n)-l) (12) 

where 

It is s easy to veriiy that a{n) satisfy the recurrence equation 

from which it can be proved by induction that a(n) = n/ D. Replacing into Eqs.(6) and (12) we finally obtain 



where in the last step we have used that 



- (-ir^ = -2 (15) 

n=l 



„=D/2 
n=l 



IV. FINITE TEMPERATURE EQUILIBRIUM PROPERTIES 

We now present some numerical simulation results at finite temperature. Monte Carlo simulations were performed 
using heat bath dynamics. Wc first analyze the frustrated case _D = 4L / = 1, 2, . . .. In Fig. 2 we show the results 
for the internal energy per spin for system sizes _D = 4, 8, 12. The different curves were obtained by starting at high 
temperature with random or ferromagnetic initial conditions (there were no difference between both cases in the final 
results) and slowly cooling, that is, the initial configuration at every subsequent temperature was taken as the last 
configuration of the previous temperature. At every temperature we left the system thermalize during a varying 
transient period, which was higher at low temperatures. 

For low values of D we see that the energy converges to the value —0.5. Direct inspection of the final configurations 
obtained at T = shows always the conjectured state composed by two antiferromagnetic domains. Notice from 
Fig. 2 that there seems to be a change of behaviour around T ^ 0.2. An analysis of the energy distribution at low 
temperatures showed that this effect is related to a gap between the ground state and the first excited states that 
decreases as D increases and hence, it is a finite size effect. Both the magnetization and the staggered magnetization 
are always 0{N~^/'^) at any temperature and also the specific heat does not show any type of anomaly. So the system 
does not (globally) order at any temperature. 
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FIG. 2. Monte Carlo calculation of the mean energy per spin vs temperature for an antiferromagnetic Hamiltonian with 
D = 4,8, 12. 

We next performed Monte Carlo simulations for the antiferromagnetic system when D/2 is odd. 

In figure 3 we show the absolute value of the staggered magnetization, the energy per spin and the staggered suscep- 
tibility vs temperature for D = 6, 10, 14. These results not only confirms that the ground state is antiferromagnetic, 
but also show that the system undergoes a second order phase transition at a finite critical temperature Tc « 1. The 
comparison of the magnetization curves with the solution of the Curie- Weiss equation m = tanh{l3m) suggests that 
mean field behaviour is exact in this case at any temperature. 
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FIG. 3. Monte Carlo calculation of the mean staggered magnetization per spin, energy per spin and staggered susceptibility 
vs temperature for an antiferromagnetic Hamiltonian with D = 6, 10, 14. 



V. THE FERROMAGNETIC MODEL 

Finally we consider the ferromagnetic version of the model, that is, the model defined by the Hamiltonian (1) with 
J(r) = — jv^- This system is non-extensive, in the sense that the quantity 

diverges for N ^ oo, and therefore the thermodynamic limit is not defined [18,19]. In this situation the magnetization 
per spin m{T) and the internal energy per spin u{T) are expected to scale as [20,21] 

m(T) ~ m'{T/4>) 

u{T) ~ cl>u'{T/<i>) (18) 

when TV ^ 1, m'{x) and u'{x) being scaling fimctions independent of N. In the case of d-dimensional hypcrcubic 
lattices with interactions that decay with the distance r between spins as 1/r" the quantity is proportional to [21] 

N* = ^— (n^-°'/'^ - l) (19) 

1 — a/a V / 
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This scaling factor generalizes the traditional scaling J J/N of the couplings in the Curie- Weiss model, i.e., a fully 
connected Ising spin system with distance-independent interactions, which corresponds to the a = case. For a > d 
the scaling factor (19) becomes independent of in the thermodynamic limit and the system is extensive, while for 
a ^ d (this corresponds, for instance, to the case of three dimensional systems wtih dipolar interactions that decay 
as l/r^ for r ^ 1) we have that N* \n.N. Thus, the behavior of the present hypercell model is expected to be 
analog to that of the borderline case between short and long range interactions a = d va hypcrcubic Bravais lattices. 

The data collapse of the rescaled quantities for different values of D shown in Fig. (4) verify the scaling (18). More- 
over, the comparison of the magnetization curves with the solution of the Curie- Weiss equation m = tanh{m/T*) 
(where T* = T / (j)) shows that mean field behavior is asymptotically exact for a system of this type when the thermody- 
namic variables are properly rescaled, as was previously conjectured for general systems with long-range ferromagnetic 
interactions [22]. Also the specific heat calculation shown in Fig. (5) displays a convergence to a discontinuous curve 
at the rescaled critical temperature when N ^ oo, consistently with the mean field behaviour. 
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FIG. 4. Monte Carlo calculation of the mean magnetization per spin and energy per spin vs rescaled temperature T* for a 
ferromagnetic Hamiltonian with D = 4,6, 8, 12. 
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FIG. 5. Monte Carlo calculation of the specific heat vs rescaled temperature T* for D = 4,6, 8, 10. 
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VI. DISCUSSION 



We analyzed the statistical mechanics of a long-range model defined on a D-dimensional hypercube, which is related 

to the complexity of Boolean functions. Besides its interest in Boolean functions complexity and neural networks fields, 
the model exhibits several features that contribute to the understanding of general long-range systems. Thus, the 
statistical mechanics of the model is of interest on its own. 

In the antiferromagnetic (AF) version of the model, the particular form of the coupling terms in Eq.(2) actually 
define two different models, according to the value of D. When D/2 is an odd number the model presents a low 
temperature AF phase and undergoes a mean-field like second order phase transition at a finite critical temperature 

= 1. The most interesting situation appears when D/2 is and even number. In this case long range AF order is 
suppressed at any finite temperature, while the ground state is composed by two equal size anti-aligned AF domains; 
this state becomes infinitely degenerated in the thermodynamic limit, but its entropy per spin goes to zero. This result 
shows that the frustration produced by the interplay between long range interactions and lattice topology can produce 
large ground state degeneracies without the generation of a "rough" free energy landscape. This result is of importance 
in spin glasses studies, where models defined on an hypercube are \iscd to analize the influence of large dimensionalities 
[11,13]. We have shown that this type of lattice can introduce effects of large ground state degeneracies which are 
independent of the presence of disorder, at least in the case of long range interactions. Although the contribution of 
such effects to the entropy per spin goes to zero in the thermodynamic limit, its importance could not be negligible for 
finite systems, as in the case of numerical simulation studies. It is worth to stress that the antiferromagnetic model 
do not have a well defined thermodynamic limit for general dimension of the hypercell D. The system exhibit a well 
defined thermodynamic limit only when the dimension D is restricted to a particular set of values, that is, to be a 
multiple of four or not. This peculiar behavior illustrates the type of subtleties that may arise when working with 
long range frustrating interactions in hypercubic cells. 

Finally, in the ferromagnetic version of the model we verified a previous conjecture [22] for non-extensive (ferromag- 
netic) systems, that is, all the thermodynamic behaviour of the model can be exactly described at any temperature 
by mean field theory if the thermodynamic variables are properly rescaled with the system size. 
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